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Abstract 

We show how the generalized Lambert series ^2 n>1 xq"/(l — xq n ) can be 
computed with Theta convergence. This allows the computation of the 
sum of the inverse Fibonacci numbers without splitting the sum into even 
and odd part. The method is a special case of an expression for the more 
general series ~^2 n>0 t n / (1 — xq n ), which can be obtained from either the 
Rogers-Fine identity or an identity by Osier and Hassen. 

1 The identity 

Fine's function F(a,b,t) can be given |5[ p.l] as 

(l~aq)(l-aq 2 )(l~aq 3 ) ■■■ (l-ac^l „ 



(1 - bq) (1 - bq 2 ) (1 - bq 3 ) ■ ■ • (1 -bq n ) 

Denning (a; q)o = 1 and (a; q) n = (1 — a) (1 — a q) (1 — a q 2 ) ■ ■ ■ (1 — ag n_1 ) (for 
n G N + ) we can write 

(a q)n , n 



F(a,b;t) = J2 



n>0 



(bq; q)r 



The Rogers-Fine identity ||, relation (14.1), p. 15] is 

(i - 1) F( a , b-,t) = y ( a q : q\(atq/b-q) n q2n+1 bn tn . 

frj; (bq;q)n(tq;q)n 



Replacing b by a q gives, for the left side, 



(i-0£fW*» = (i-t)£ T 
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For the right side we get 

y (a g; g)n (*? g)n (1 _ a t g 2n+U {aq) n f n = 
^ {aq 2 ;q) n {tq;q) n 

Now we divide both sides by (1 — a g) (1 — t) to obtain 

f n i 

y —^T = Y 7 TTw r {I -at q 2n+1 ) (aq) n t n q n 

^l_ ag n+i Z-^ (1 - a (1-t q n ) nq > 

Finally replacing a by x/q gives the desired identity 

^1-xf 4^ (1 - a;<7") (1 ~ tq n ) w 

Both sides are convergent for \q\ < 1, |x| < 1, and |t| < 1. The right side reveals 
symmetry in x and i which is also observed as relation (4.3) in Jl3| , p. 655]: 

j-n n 

y — — = y — — (2) 

1 - x q n ^ 1 - t q n 

n>0 y n>0 ^ 

Setting x = q, t = q in (Q), and multiplying by q gives the Clausen's g p.95] 

n>l y n>l y 

We denote the rate of convergence of the series on the right side as Theta 
convergence. (This is faster than linear but not what is usually called super- 
linear convergence, which should really be called exponential convergence.) 



2 Summation of inverse Fibonacci numbers 

The computation of ^ n>1 l/F„ 3.359885... where F = 0, F\ = 1, and 
F n = F n -i + F n -2 for n > 2 (the Fibonacci sequence) has been considered 
in various publications (see, for example, Q, S, ||, and |J). The following 
identities for the sums with even and odd indices have been known for a long 
time. Let a = 1/2 (1 + y5), /3 = — 1/a = 1/2 (1 — \/S), then (even indices, see 
@ P-468]) 

V -J- = y/E (L{f3 2 ) - L(/3 4 )) « 1.5353705 .. . 
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and (odd indices, see pp.94] 

V5 



^-^4W^ 9 ^) « 1.8245151. 



n>l 



where 63(g) = 1 + 2 J^ n>1 e? n • Both relations allow computation with Theta 
convergence. Other expressions for the sums are 

n{n+l) 



and 



^ 1 _^P n ^ n+i > 

n>l n>l 
n>l \n>0 / 



The search for a relation that does not require the splitting between even and 
odd indices (and works for more than just the Fibonacci numbers) has been 
the main motivation for this work. Gosper gave the following formula for the 
Fibonacci sequence @, p. 66]: 

y 1 =y (-1)" ( - 1)/2 (*4n+3 + (-!)" F 2n+ 2) (4) 
^ — ' F n -Pan+l ^2n+2 Gi G3 G5 • • • G2n+1 



n>l 



where G„ = 2F„_i + F„ (Lucas numbers). Note the last term in the denomi- 
nator has been corrected from Gin-\ to Gm+\- 

We will need a specialization of (|l]), often called generalized Lambert series. 
First set t = q and multiply by x to obtain 

xq n = y, (l-xq 2n+1 ) ^ 2 +n 



Replace x by x q to obtain the desired form (note the new limits of summation) 

q n ^ {l-xq 2n ) 



v Z^i_ xq n (1 -xq n )(l-q n ) H W 



The computation of J2 n >i V/« where /o > 0, fi > 1, and f n = mi f n -i + 
m% fn-2 for n > 2 has been treated by Horadam We restrict our attention 
to cases where fo = 0, fx = 1 (Fibonacci type sequences), mi > 1, m-z 7^ 0, 
and A = m\ + 4 7712 > 0. The special case mi = 7712 = 1 corresponds to the 
Fibonacci numbers. The roots of the characteristic polynomial x 2 — (mi x + m 2 ) 
are a and j3 where 

a= 1/2 (mi + Va) , P = 1/2 (mi-VA) 
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The following identity, given by Horadam || relation (4.17), p. 108], together 
with relation (||), is the key to fast computation of J2 n>1 1//«: 



n>l 



3 A few more relations 

We mention several identities observed during this work. 

3.1 Alternate relation for fast computation 

Fine's second transformation is || relation (12.2), p. 13] 

y (aq^ n = y {b/q;q) n ( „ 2+ „ )/2 
k ' ^ o (bq;q) n ^{bq;q) n {tq;q) n K 

Setting b — aq, then replacing a by a/q and dividing by (1 — a) (1 — t), and 
finally replacing a by x gives an alternate identity for fast computation: 

V l l -V (g ;< ?)" ( Tf \n „(n 2 -n)/2 m 

The symmetry between a; and t is evident from the right side, giving (^|) once 
more. 

3.2 Another derivation of the main identity 

Identity (Q) can alternatively be deduced from the (remarkable!) relations (6.5), 
(6.6), and (6.7) by Osier and Hassen Jlj| p. 659], namely 

~,n ~d (a n+b) on „fc (c n+d) 

a q = p q = 

Z—/ \ _ R a c(a n+b) Z^i \ — nn a ( c n+d) ^ ' 

n>0 H y n>0 H 

( oH 1 (3 n q( an + b ) ( cn + d ) a n+1 p* n g( a (™ +1 )+ b ) ( cn + d ) \ 

= ( I _ p q c (a n+b) 1 1 _ aq a(cn+d) ) ,S "' 



n>0 



by setting a = c = 1, = d = 0, and replacing a by x and /3 by t. 
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3.3 The symmetric relation (0) from an identity by Fine 

Fine's relation (16.3) §, p.18] is 

cv h +\ i a 1il)oo {b/a;q) n (aq) n 
F(a,b,t)= ^ 



( b ^l)o° n> (9;q)n i-tq n 

The following can be obtained by setting b = aq, dividing by (1 — aq), and 
finally replacing a by a;: 

^ 1 - X q n+1 ^ 1 - t q n 

n>0 n>0 

Replacing x by x/q we get relation (^) once again. 

3.4 More symmetric relations 

An identity resembling ^ is 

V - = V 



^ o (.t; 9 )„ +1 ^ Q (t;q)n+l 

This follows from an identity given by Gosper || 

(t; q)oo(x; q)oo \^ (t;q) n (x;q) 

(?;?)« ^ (x;q)„ +1 ^ (q;g)„ 

and noting that both the right side and the factor in front of the sum on the 
left side are symmetric in x and t. 



Setting A = B = C = D = lin relations (|8aJ) and (8b) (and again replacing a 
by x and j3 by t) we obtain 

n n 4-n n 

4^ 1 - * <7" ^ 1 



n>l n>l y 



v-^ X n t n (J"' v-^ X n t n 
> h X > 

^ 1 - 1 q n ^ l-xq n 

n>l y ri>l y 

x™ t n q n ' 2 ^ x n t n o n C" +1 ) 

> — +0 " 

^ l-xq n ^ l-tq n 
n>l y n>l y 



^ (1 - x o") (1 - f a") 



n 4- n Jn 



^(l-xq")(l-tq n ) 
The last form reveals the symmetry in x and t. 
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3.5 Another derivation of the relation for the generalized 
Lambert series 



Relation (g) can be obtained from the following relation given by Knuth (at- 
tributed to J. R. Wrench, Jr.) p. 644, solution to exercise 5.2.3-27]: 



E a n Q _ \ ^ 
\ _ nil 



a n + E ("« + a n+k) 



n>l 1 n>\ 

In the following (well-known) relation 

x q" 



k>l 



Ex q \ ^ -h 

l-xq n ~ ^ 1 - q ri 

n>0 n>l 



replace x by g x to get 



1 _ xq n 1 _ «n 



TI>1 



Now set a„ = x™ to obtain 



= E 



n>l 



l + ]T(l + x fe ), 



fe>i 



x n q n 



which gives identity (|j) in the form 



L(x,q) = J2 



n>l 



E 

n>l 



X q 



fe>0 



fc>l 



x n g" 



1 - q n 1 - x q r ' 



x n q n 



3.6 Gosper's matrix-product for Clausen's identity 



Gosper defines || matrices K(k,n) and N(k,n) as 

r ,n+2k+l 

I 



K(k,n) 



q ..T~~^ q (l_ q 2k+ny((i _ g fe) (x _ 



<? fc g/(l - q fe+ ™)' 
1 



N(k,m) 

These matrices satisfy 

N(k, n) ■ K(k, n + 1) = K(k, n) ■ N(k + 1, n) 

Writing K(k, oo) for lim n _ i . 00 K(k, n) and N(oo,n) for limfc_>. 00 iV(fc, n) we ob- 
tain (^) as upper right entries on both sides of 

[J N(l, n) ■ 11 K(k, oo) = n 0) • J] tf(oo, n) 



ra>0 



fe>l 



k>l 



n>0 



G 



3.7 Jordan's identity for the bilateral series 

For the bilateral series (Jordan-Kronecker Function), defined as 

t n 



x q" 

n— — oo 



we have f(x, t) = f(t, x) and the following two relations given in @, chap.3] 

/(*,*)= E ^^^(+1 + ^- + ^-) 

v ; ^ V l-xq" l-tq n 

n— — oo 



+ 00 , 

f(x,t) = V g" 2 o; n i"(-l+ + 

^ v ' y ^ y V 1 - x a" l-ta n 



n— — 00 



As pointed out by Shaun Cooper, the latter is already given in |10[ rel.(22), 
p. 450]. Both forms are equivalent to 



+ °° C\ - t-rn 2n \ 

f( X ,t)= y 7 [ , q ' , x n ef 
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